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Abstract 

The purpose of this paper is to construct a generalized r-matrix structure 
of finite dimensional systems and an approach to obtain the algebro-geometric 
solutions of integrable nonlinear evolution equations (NLEEs). Our starting 
point is a generalized Lax matrix instead of usual Lax pair. The generalized 
r-matrix structure and Hamiltonian functions are presented on the basis of 
fundamental Poisson bracket. It can be clearly seen that various nonlinear 
constrained (c-) and restricted (r-) systems, such as the c-AKNS, c-MKdV, 
c-Toda, r-Toda, c-Levi, etc, are derived from the reduction of this structure. 
All these nonlinear systems have r-matrices, and are completely integrable 
in Liouville's sense. Furthermore, our generalized structure is developed to 
become an approach to obtain the algebro-geometric solutions of integrable 
NLEEs. Finally, the two typical examples are considered to illustrate this ap- 
proach: the infinite or periodic Toda lattice equation and the AKNS equation 
with the condition of decay at infinity or periodic boundary. 

Keywords Lax matrix, r-matrix structure, integrable system, algebro-geometric 
solution. 
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1 Introduction 



Completely integrable systems are widespreadly applied in field theory, fluid mechan- 
ics, nonlinear optics and other flelds of nonlinear sciences. The new development of 
integrability theory can be roughly divided into three stages. The flrst one was the 



direct use of Lax equations for some systems such as the Calogero-Moser system |33 



and the Euler rigid equation , which were allowed for integration. These systems 
were not amenable to the classical techniques of integrating Hamiltonian equations 
of motion. The second one was the so-called 'algebraization', i.e., the tools of Lie 
algebras, Kac-Moody algebras were used to sysmetically construct a large class of 
soliton equations and integrable systems ^ 0, and simultaneously present the 
Lax representations of soliton equations, Hamiltonian structures. The third one is 
being developed and witnessed through the use of nonlinearization technique [§ to 
generate finite-dimensional integrable systems. These systems can be the Bargmann 
system, the C. Neumann system [|1^, the higher-order constrained flows or sym- 



metric constrained flows 0, and the stationary flows of soliton equations |53 



Indeed, with the help of this method, many new completely integrable systems were 
successively found [p!2| |35| . In this way, each integrable system is generated through 
making nonlinearized procedure for a concrete spectral problem or Lax pair, and it 
has its own characteristic property. Then a natural question is whether or not there 
is a unifled structure such that it can contain those concrete integrable systems? 
Recently, the study of r-matrix for nonlinear integrable systems brings a great hope 
to solving this problem. 



Semenov-Tian-Shansky ever gave the deflnition of r-matrix |47], and used the 



r-matrix to construct Lie algebra and new Poisson bracket |^ in a given Lie algebra 
and corresponding coadjoint orbit. The main idea of Semenov-Tian-Shansky and 
Reyman was how to obtain the new Poisson bracket from a given r-matrix and an 
element of Lie algebra. Here our thought is how to present r-matrix structure from 
a given Lax matrix and the standard Poisson bracket: 

L, {■,-}=^ r — matrix. 

In the present paper, we give a sure answer for the above question. We propose 
an approach to generate flnite dimensional integrable systems by beginning with 
the so-called generalized Lax matrix instead of usual Lax pair. Another main result 
of this paper is to deal with the algebro-geometric solutions of integrable nonlinear 
evolution equations (NLEEs). It is well-known that the ideal aim for nonlinear 
equations is to obtain their explicit solutions. According to the nonlinearization 
method, solutions of integrable NLEEs can have the parametric representations 



ro or involutive representations and also have numeric representations in the 
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discrete case |4||. However, these representations of solutions are not given in an 
explicit form. Thus, an open question is how to obtain their explicit forms. In the 
paper we would like to give solutions of integrable NLEEs in the form of algebro- 
geometric G-functions. 

The algebro-geometric solutions for some soliton equations with the periodic 
boundary value problems were known since the works of Lax |2^, Dubrovin, Mateev 
and Novikov [^. Similar results for the periodic Toda case were obtained slightly 
later by Date and Tanaka fl^. Afterwards, the relations between commutative 
rings and ordinary linear periodic differential operators and between algebraic curves 
and nonlinear periodic difference equations were discussed by Krichever The 
technique they adopted is the Bloch eigenfunctions, the spectral theory of linear 
periodic operators, and some analysis tools on Riemann surfaces. 

In the first example of this paper, we give the algebro-geometric solutions of the 
Toda lattice equation in the infinite or periodic case. Our method is a constraint 
approach connecting finite dimensional integrable systems with integrable NLEEs 
instead of usual spectral techniques and Bloch eigenfunctions which are often avail- 
able to the periodic boundary problems. The results with the periodic boundary 
conditions are included in ours. In the second example, we consider the well-known 
AKNS equation. The Ablowitz-Kaup-Newell-Segur (AKNS) equations are a very 
important hierarchy of NLEEs in soliton theory. It can turn out that the KdV, 
MKdV, NLS, sine-Gordon, sinh-Gordon equations etc. All these equations are solv- 
able by the inverse scattering transform (1ST) |2^, and usually have A^-soliton solu- 
tions 1^ . But the algebro-geometric solutions of the AKNS equations have not been 
obtained since then. We shall solve this problem by using our constraint procedure. 
The considered AKNS equation is under the case of decay at infinity or periodic 
boundary condition. 

The whole paper is organized as follows. We first introduce a generalized Lax 
matrix in the next section, then construct a generalized r-matrix structure and a 
generalized set of involutive Hamiltonian functions in section 3. All those Hamilto- 
nian systems have Lax matrices, r-matrices, and are therefore completely integrable 
in Liouville's sense. In section 4 it can be seen that various constrained (c-) and 
restricted (r-) integrable flows, such as the c-AKNS, c-MKdV, c-Toda, r-Toda, c- 
WKI, c-Levi, etc, can be derived from the reductions of this structure. Moreover, 
the following interesting facts are given in sections 5, 6, 7, respectively: 

- Several pairs of different integrable systems share the same r-matrices with 
the good property of being non-dynamical (i.e. constant). In particular, a discrete 
and a continuous dynamical system possess the common Lax matrix, r-matrix, and 
even completely same involutive set. Additionally, on a symplectic submanifold 
integrability of the restricted Hamiltonian flow (for continuous case) and symplectic 
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map (for discrete case) are described by introducing the Dirac-Poisson bracket. They 
also have the same r-matrix but being dynamical. 

- A pair of constrained integrable systems, produced by two gauge equivalent 
spectral problems, posesses different r-matrices being non-dynamical. 

- New integrable systems are generated through choosing new r-matrices from 
our structure, and the associated spectral problems are also new. 

In the last section, as a development of the generalized structure, through considering 
the relation lifting finite dimensional system to infinite dimensional system and 
using the algebro-geometric tools we present an approach for obtaining the algebro- 
geometric solution of integrable NLEEs. To illustrate the procedure we take the 
periodic or infinite Toda lattice equation and the AKNS equation with the condition 
of decay at infinity or periodic boundary as the examples. 

Before displaying our main results, let us first give some necessary notation: 
dp A dq stands for the standard symplectic structure in Euclidean space R^^ = 
{{p, q)\p = (pi, . . . ,Pn), <? = (<?!, • • • , Qn)}] < ■, ■ > is the standard inner product in 
R^; in the symplectic space (-R^^, dp A dq) the Poisson bracket of two Hamiltonian 
functions F, G is defined by 

^1 _ Vr— — - — — ) - dF^ dG 

i^i 9qi dpi dpi dqi dq ^ dp dp' dq ^ 

1 and (8> stand for the 2x2 unit matrix and the tensor product of matrix, respectively; 
Ai, Atv are arbitrarily given distinct constants; A, /i are the two different spectral 
parameters; A = diag{Xi, Xjy), Iq =< q,q >, Jq =< p,q >,Kq =< p,p >, 
Ji =< Ap,p >< Aq,q >, Ji =< Ap,q >, oq, ai = const.. Denote all infinitely 
times differentiable functions on real field R by C°°(i?). 

2 A generalized Lax matrix 

Consider the following matrix (called Lax matrix) 

A{X) B{ 
C(A) -A{\) 



m-iT. (2.1) 



where 



N 

A{X) = a_2(/i,Ji)A-2 + a_i(Jo)A-^ + ao + aiA + E-^, (2.2) 

BiX) = 6_i(/o,Jo)A-i + 6o(Jo)-Et^' (2.3) 

j=i ^ - 

N 2 

C(A) = c„i(Jo,iro)A-i + co(Jo) + ET^- (2-4) 
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with some undetermined functions a_2! d-i, b^i, c_i, Bq, cq E C'^{R). 

Now, in order to produce finite dimensional integrable systems directly from the 
Lax matrix (2.1), we need an inevitable Assumption (A): {A(A), A(/i)}, 
{A(A),5(/i)},{A(A),C7(/i)},{B(A),B(/i)},{fi(A),C(/i)}, and {C(A),C(/i)} are ex- 
pressed as some linear combinations of A{\), v4(/x), -B(A), -B(/i), C(A), C(/x) with the 
cofficients in C°°{R). Then we have 

Lemma 2.1 C/nc?er Assumption (A), L{X) only contains the following cases: 

1. If a-2 7^ const., ao — bo — Co — ai — 0, a_i = — Jo, = -^o, andc-i — —Ko, 
then a_2 satisfies the relation Ii — (Ji + a_2)^ + /(a_2); a-2 = const. ^ anc? 
do — bo = Cq = Qi = 0, ^/len a_i = —Jo, ^-i = -^Oi c_i = -Kq, or a_i = 
const., 6_i = /q + fi{Jo), c_i = — i^^o + 9'i('^); where fi, gi satisfy the relation 
fi9i = -Jq - 2a_iJo + const.. 

2. a„2 = fl-i = b^i = c_i = bo = Cq = ai = 0, and ao = const. 

3. a_2 = 6-1 = ao = bo = Co = ai = 0, c_i = —Kq, and a_i satisfies ^^^^ 7^ 0. 
^. a_2 = a_i = 6-1 = c_i = 60 = 0.1 = 0, ao = const., and co 7^ 0. 

5. a_2 = c_i = 60 = ai = 0, a_i, ao = const., b-i = lo + g{Jo), and co satisfies 
^M(^o)) = -2ao. 

^. a_2 = c_i = Oo = &o = Co = Oi = 0, a_i = Jo + const., and 6_i = /o- 

7. // a_2 = ao = 60 = Co = ai = 0, ^/len ^/lere are the following five subcases: 

(7.1) a_i = const., c_i = -Ko + /2(Jo), and 6_i = /q +52(^0); 

(7.2) a_i = - Jo, 6_i = /o, and c_i = Kq; 

(7.3) o_i = — Jo + const., and 6_i = 6_i(Jo),c_i = c_i(Jo) satisfy ^(6_ic_i) = 

2a_i; 

(7.4) a_i = — Jo + const., 6_i = /o, and c_i = c_i(Jo); 

(7.5) a_i = — Jo + const., c_i = —Ko, and 6_i = 6_i(Jo). 

8. a_2 = a_i = 6_i = c_i = 0, ao, ai = const., bo 7^ 0, Co 7^ 0, anc? 60, Cq satisfy 
the relation -^{boCo) — — 2ai. 

9. a_2 = a_i = 6-1 = c_i = 0, co, ai, ao = const., and bo 7^ 0. 

10. a_2 = = Co = ai = 0, a_i, ao = const., c_i = —Ko + h{Jo), and bo 
satisfies the relation ^(^o^(<^)) = — 2ao. 

The above all functions f, g, h, fi, g^ {i = 1,2) are in C°°{R). 

Proof Through some calculations we have 

{A{X),A{i,)} = 2^<Ap,p>{^BiX)-^B{^,)) 
+2^ <Aq,q> (Ac(A) - i^C(/x)) 

+2-^(^ - ^)(< Ap,P > (6_i - /o)- < Ag,(Z > (c_i + i^o)) 
- ^ < ^P^P > ^0+ < Ag,g > Co), 

ail A^ 
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{B{X),B{f,)} = 2^(li?(A)-^5(M)) + 2^(S(A)-5(/x)) 

1 1 db-i dbo dh-i dbo 

{C7(A),C(/.)} = 2^(-1c(A) + ^C(m)) + 2^(-C7(A) + C(/x)) 

{A(A),B(/.)} = ^(-B(m) + B(A))-|^S(m) 

2 db-i 2n Ba-2 , 5a_2 . ^ 

- + ^ «r < ^"'^ > 

-2(2 < Ag, g > ^(Jo + a_i) + - ?'-i))A-2 

-2 -60^ + 2ao < Ag, g > ^ 

, da_i , 1 , 1 , . da-2 , 9 9 

+26o-TyiA-i - 26o^/x-^ - 4ai < Ag,g > — -iA-^', 
aJo a/o oil 

2dc-i^,^. 211. da-2 ^, . „9a_2 . .... 

+2 < Ap.p > (2|^^ J, + - 2a_,55_l)A-V-' 

aiio oil oKq oJi oil 



,dc-ida^i^^ dc-i da-i _i 

oKq dJo oKq dJo 



-2(-2 <Ap,p> ^(Jo + a-i) + ^{Ko + c_i))A-2 

-2 CO— ^ + 2ao <Ap,p> ^ A'^^u 

-2co— -A -2co^— /u ^ -Aai<Ap,p> ——X V, 
aJo aiio c/Ji 



4 4 96 1 

{B(A),C(m)} = ^(-A(m) + ^(A)) + -^^(m) 



fi oKq fi oJq dJo X oJq dJo 

+ l)«-2/^-'A-^ + 4(^ + l)a_2A- 
olo oKo 

(?c_i 56-1 o^^-i ^'^-1 7 ^^-1 (^c-i 
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db^i 5c_i , db^i (9c_i , , , , 

^c_i - ^6-1 - 2^a_i + 2— ia_i + 2a_i A" V"^ 
oJo oJo oIq oKq 



,db() dcn , , ,c^c_i 1 , 1 

-2 -^co + -f 6o + 2ai + 4— iai^^U - 4— ^aiA-V- 
aJo aJo c'.f^o cf/o 

According to Assumption (A), the terms that do not contain v4(A), v4(/i), B{X), 
B{ii), C(A), C(yu) in the above six equahties, are zero. After discussing these terms, 
we can obtain every result in Lemma 2.1. I 

3 Generalized r-matrix structure and integrable 
Hamiltonian systems 

Let Li{\) = L{X) /, L2{fi) = / ® L{fi). In the following, we search for a general 



4x4 r-matrix structure ri2(A,/i) such that the fundamental Poisson bracket ^1 



{L(A) ? = [n2(A,/i),Li(A)] - [r2i(/i,A),L2(/i)] (3.1) 

holds, where r2i(A,yu) = Pri2(A,/i)P, P = ^J2i=o'^i ® <7j) and cr-s are the standard 
Pauli matrices. For the given Lax matrix (2.1) and the Poisson bracket (1.1), we 
have the following Theorem. 

Theorem 3.1 t/n(ier Assumption (A), 

n,{X,fi) = ^—P + S (3.2) 
— A 

is an r-matrix structure satisfying (3.1), where 





/ 


2A 


da_2 2 da_i 




1 Q''-! 






dJi (A. djQ 






s = 




2A 


dJo 

<Ap,p> gj^^ 






2 86-1 
M S/o 




V 







2A 


< Ap,p > 



-2 



dbg 



'dJo 

9n-2 2 9c-l 2A fl''-2 I 2 

c?Jo c?Ji djQ 



Proof Under Assumption (A), we have 



{A{X),A{^,)} = 2^<Ap,p>(A5(A)-iiE(/.)) 



+2^<Ag,g>(Ac(A)-|^C(/i)), 
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{Bi\),Bi^^)} = 2f-i(is(A)-ii?(/i)) + 2^(S(A)-B(/i)), 
oJq jj, A dJQ 

{C{\)M^^)} = 2^(--C(A) + ^C(^)) + 2^(-C(A) + C(^)), 

9 9 dn 9 r)/) 

{A(A),5(/.)} = i-Bifi) + BiX)) - -—^Bifi) + -^B{X) 

{A(A),C(/.)} = Y^(C(/.)-C(A)) + |^C(/i) + -|^C(A) 

A — jJ, A dJo jJ OKq 

+^ ^ a/f ^'"^ ^ ^ < Ap,p > 

{5(A), = (-A(/.)+A(A)) + -^A(/i)--^A(A) 

A — jl A oIq /i aA 

which completes the proof. I 
In general, Eq. ( p. 21) is a dynamical r-matrix structure, i.e. dependent on canon- 
ical variables Pi^qi 0. 

Now, we turn to consider the determinant of L{X) 

-detL(A) = ]^TrL^{\) = A^{X) + B{\)C{X) 

= EH^X^ + Et^ (3-3) 



where 



A - Aj 



H-i = al^, (3.4) 

F_3 = 2a_2a-i, (3.5) 

H-2 = 1 + 2a-_2ao + 6_ic_i - 2a_2 < A~^p,g >, (3.6) 
= 2a_2«i + 2a_iao + 6-iCo + 6oc_i - 2a_2 < A^^p, g > 

-2a_i < A^V, 9 > -^-1 < A"V,P > +C-1 < A^^g, 17 >, (3.7) 

-f^o = al + 2a-iai + boCo + 2ai <p,q>, (3.8) 

i^i = 2ooai, (3.9) 

H2 = ai, (3.10) 
Ej = (2a_2Aj^ + 2a_iAj^ + 2ao + 2aiXj)pjqj 

+(6„iAt1 + 6o)p| - (c-iAji + co)q] - T,, (3.11) 

r,. E ,3.12, 
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Let (3.3) be multiplied by a fixed multiplier {k G Z), then it leads to 

2 k—1 N \ k IP 



where 



2 7 J ■ ^ ■ -, A — A-i 

-1 fe-1 

= X! ^l-k^^ + y^X^l-k + -^fc-l-/)A^ 
/=fc-4 /=0 

fc+2 Af \ A: p 

+ E^'-^A' + ET3f (3-13) 
^m = E^r^i' ^ = 0'1'2,... (3.14) 



which reads 

= 2a_2 < A'^-^p, q > +2a_i < A'""^^, q > +2ao < A"p, g > +2ai < A'"+^p, g > 
6-1 < X^~^p,p > +bo < K^p,p > -c_i < K^'^q,q > -cq < K^q.q > 
- E {< A''p,p >< A^q,q > - < A'p,q X A^p,q >). (3.15) 

i+j=m—l 

Because there is an r-matrix structure satisfying (3.1), 

{L'iX) « L\fx)} = [n2(A,/i),Li(A)] - [f2i(^, A),L2(/x)], (3.16) 

where 

f,,(A,/i) = ^^Lj-'=(A)Li"'(/x) -r.^^A,/.) ■Lt(A)4(/i), = 12,21. (3.17) 

A:=0 1=0 

Thus, 



4{TrL\X), TrL\fi)} = Tr{L\X) ® L\fi)} = 0. (3.18) 

So, by ( |3.13| ) we immediately obtain 

Theorem 3.2 t/n(ier assumption (A), the following equalities 

{E,, E^} = 0, {Hi, E^} = 0, {F^, E^} = 0, (3.19) 
^,j = l,2,...,iV, / = -4,...,2, m = 0,l,2,..., 

/ioW. Hence, the Hamiltonian systems {Hi) and {Em) 

r)H f)H 

{Hi) : g, = — i, p, = -^, / = -4,...,2, (3.20) 
op oq 

BE dE 
{Fra) : ^^r. = ^, Ptr. = -^, m = 0,1,2,..., (3.21) 

are completely integrable in Liouville 's sense. 

Corollary 3.1 All composition functions f{Hi,Em), f G C°^{R) are com- 
pletely integrable Hamiltonians in Liouville's sense. 
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4 Reductions 



For the various cases of Lemma 2.1, we give the corresponding reductions of r-matrix 
structure ri2(A,/i) in this section. The following numbers of title coincide with the 
ones in Lemma 2.1, i.e. the corresponding conditions are coincidental. 

Before giving our reductions, we'd like to re-stress the two "terminologies" used 
usually in the theory of integrable systems in order to avoid some confusions: one 
is ^^constrained system" , which means the finite dimensional Hamiltonian system or 
symplectic map in i?^-^ under the Bargmann-type constraint; the other ^^restricted 
system" , which means the finite dimensional Hamiltonian system or symplectic map 
on some symplectic submanifold in R"^^ under the Neumann-type constraint. In the 
future we shall follow this principle. 
1. 



ri2(A,/i) 



2A 
/i(/i - A) 



OJi /i^ oil fi^ 



S = 



( 


1 












/ 








<Aq,q> 





\ 
















: Q = 













- <Aq,q> 




















- < Ap,p> 
























1 } 




I 





< Ap,p> 











Particularly, with /(a_2) = —1, (4-V exactly reads as the r-matrix of the con- 
strained WKI (c-WKI) system / p^/ . With a_2 = const. ^ 0, (3.2) reads the r-matrix 

'^12(A,/i) 

— -P of ellipsoid geodesic flow j^W, or reads as 



n2(A,/i) 



2 2 
-P+-S, S 

jJ, — A /i 



I ^ f'l \ 

0-10 
0-100 

V -g[ ^ ) 



which is a new r-matrix structure. For simplicity, below write V'= 



2. 



n2(A,Ai) 



-P. 



(4.2) 



/i — A 

This is nothing hut the r-matrix of the well-known constrained AKNS (c-AKNS) 
system /|^. 
3. 

\ 
0-10 




n2(A,^) 



2 2 

-P + -S, S 

/i — A n 



a'„i ^ 0. 



(4.3) 



\ a'_i y 
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In particular, with a_i = —Jq Eq. (4-3) reads as the r-matrix of the constrained LZ 
(c-LZ) system 

I 

/ \ 



ri2 = rP + c'^S, S 

/i — A 



10 

0-10 



(4.4) 



With Co = —2y/J^, Eq. (4-4) reads as the r-matrix of the constrained Hu (c-H) 
system / |7^/ . 
5. 

( 

n2(A,/i) 



-P + 5, S 



-g' Q \ 
c'o 

-V 



(4.5) 



/i — A 

V y 

VFzt/i 6_i = /q, Co = const., Eq. (4-5) reads as the r-matrix of the constrained Qiao 
(c-Q) system jg^/ . 
6. 

/ 1 \ 



ri2(A,/i) 



r/izs zs a new r-matrix. 
7. 

(7.1) 



2 2 

-P + -5, 5 



— A /i 




10 
1 



(4.6) 



'^12(A,/i) 



2 2 

-P + -S, S 



^ — X fi 



f g'^ 0\ 

/2 -1 
0-100 





(4.7) 



This is also a new r-matrix. 
(7.2) 



n2(A,/i) 



2A 



(4.8) 



/i(/i-A) 

This is the r-matrix of the constrained Heisenberg spin chain (c-HSC) system j^. 
(7.3) 

/ 1 \ 




2 2 
ri2(A, /i) = -P 5, S* 

11 — A yU 




\ c'_i I j 



(4.9) 
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With b'_i = —1, c'_i = 1, Eq. (4-9) becomes the r-matrix of the constrained Levi 
(c-L) system 

^4) 



ri2(A,/i) 

This is a new r-matrix. 
(7.5) 

n2(A,/i) 



2 2 

-P - -S, S 

/X — A n 



( 


1 











\ 























1 




















1 


/ 



2 2 

-P - -S, S 

ji — X /i 



/ 


1 











\ 










1 






















V 











1 


/ 



(4.10) 



(4.11) 



This is also a new r-matrix. 



H — X 



P + S, s 



/ 





26'o 








\ 

































V 














/ 



(4.12) 



With bo = Cq = \/Jo, ai = — |, Eq. (4-12) reads as the r-matrix of the constrained 
Tu (c-T) system [1^] . 
9. 

/ 1 \ 

2 

ri2(A,/x) 



fi ~ X 



-P + S, S = b' 




-1 




(4.13) 



With bo = — Jo, Eq. (4-13) reads as the common r-matrix of the constrained Toda 
(c-Toda) system (a discrete system) and the constrained CKdV (c-CKdV) system (a 
continuous system), which can be seen in subsection 5.1. 
10. 

b'o \ 

2 ^h' 0-^0 

n2(A,/i) = - -P + S, S 



— X 



V 



^h' 

-6^ 
-i/i' 



(4.14) 



With h = const., bo = 0, Eq. (4-14) reads as the r-matrix of the constrained MKdV 
(c-MKdV) system ||f . 

Proof For simplicity, we only present the proof in cases 1 and 10, other cases 
are similar. 
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Case 1: With a_2 7^ const. ^ the matrix S becomes 



V 



2A do--2 _ 2 
"Jj^ dJi iJ. 



-^<Ap,p> * 



ail 



^ <Ap,p> ^ 



_2 






\ 



2A 8°-2 _ 2_ 



Substituting S into (3.2) and sorting it, we can obtain (4.1), where a_2 satisfies the 



relation Ji = ( Ji + a_ 



/(a_2)) for any /(a-2) £ C°°(i?). Particularly, choosing 



/(a-2) = -1 yields a_2 = < Ap,p >< Aq, q >- < Ap,q > . Thus Eq. (4.1) 
reads 



1 



— A) /^^ jj,'^ < Ap,p >< Aq, q > 
while the corresponding Lax matrix L{X) becomes 

hi <q,q>X-^ 



Q, (r - WKI) 



< p,p > X ^ 



-hi J PiX- Xj V p'j -Pjqj 



where 



^11 = (\/l+ < >< Ag, g >- < Ap, g >)A ^- < p,q > X \ 
Set an auxiliary matrix Mi as follows 



Ml = Ml (A) = 



A 

<Ap,p> 



<Aq,q> \ 
^/l+<Ap,p><Aq,q> 

-X X 



^l+<Ap,p><Aq,q>' 

then the Lax equation 

L, = [Mi,L] 

is equivalent to the following finite dimensional Hamilton system {^/H^i): 

djH-4, 



H-^) : 



<Ap,p> A ^ _ 



with 



^i?-4 = a_2 = - < Ap, g > +^^1+ < Ap,p >< Ag, g >, 



which is obviously integrable by Theorem 3.2. 
Let 



< A/;, p > 



V'l+ < Ap,p >< Aq,q>' 



< Ap,p >< Aq,q>' 



(4.15) 



(4.16) 
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then is nothing but the WKI spectral problem [51 

/ -A 



Vx 



Xv 




with the above two constraints ( |4.1(j| ), A = Xj, and y = {qj,pj)'^, j = 1, . . . , N. That 
means that (r — WKI) is the r-matrix of the integrable constrained WKI (c-WKI) 
system (|4.15| ). 

Other subcases in case 1 can be similarly proven. 

The proof of Case 10 can be found in ref. |3^. In this case, the corresponding 
constrained system is reduced to the well-known MKdV spectral problem [Q. I 

Remark 4.1 From the above formulas (4.1)-(4.14), the r-matrices of cases 
2, 6, and (7.2) are non-dynamical. But in fact, for other cases we can also obtain 
non-dynamical r-matrices if choosing some special functions, for instance, in Eq. 
(4.3) setting a_i such that a'_i = const, leads to a non-dynamical one. Of course, 
we can also get dynamical r-matrix, for instance, in Eq. (4.4) choosing cq = —2y/jQ 
yields a dynamical one. 

Remark 4.2 Equalities (4.1)-(4.14) cover most r-matrices of 2 x 2 constrained 
systems. But among them there are also some new r-matrices and finite dimensional 
integrable systems like cases 6, 7.1, 7.4, and 7.5 (also see section 7). Their Lax 
matrices are altogether unified in Eq. (2.1). So, quite a large number of finite 
dimensional integrable systems are classified or reduced from the viewpoint of Lax 
matrix and r-matrix structure. 



5 Different systems sharing the same r-matrices 

In the above r-matrices, we find some pairs of different integrable systems sharing 
the common r-matrices. Now, we present these results as follows. 



5.1 The constrained Toda and CKdV flows 

Let us consider the following 2x2 traceless Lax matrix (corresponding to case 
9 in section 4) 

tTC _ TTCf^ ^ _ I <p,q>^ 



( Atc{X) Btc{X) 
V CtcW -ArciX) 

where 



(5.1) 



Lo = Lo(A,p,g) = f:-^f ). (5.2) 
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The determinant of Eq. leads to 



N j^TC 



-X^+ <p,q> X + 2Hc + J2 



j=i ^-^j 



(5.3) 



= ^jPj^J - P'j- <P^Q>9j- J = 1, 2, N, (5.4) 
where Tj is defined by ( |3.12| ) and the Hamiltonian function Hq is 



Hc = <p,p> +^ <Aq,p> <q,q><p,q> . 



(5.5) 



Viewing the variables q and p as the functions of continuous variables x, then we 
have the following Hamiltonian canonical equation (He)'- 



Px 



8Hc 
dq 



-\Kp + 1 <q,q>p+ <p,q> q, 



dHc 



Qx = = -P + l^q- I < q,q> 

which is nothing but the coupled KdV (CKdV) spectral problem 



(5.6) 



-1 



V 



with the two constraints (Bargmann-type) 

u=<q,q>, v=<p,q>, 



(5.7) 



(5.^ 



A = \j and ip = {Pj, qjY ■ So, {He) coincides with the constrained CKdV (c-CKdV) 
fiow. 

Let us consider endowing with an auxiliary 2x2 matrix Mt as follows 



9 

_1 \-<q,q> I ' 

a a 



=< Aq,q> - <p,q> - < q,q>'' 



(5.9) 



Then, we have the following theorem. 

Theorem 5.1.2 The discrete Lax equation 



TC 



'TC\i 



(5.10) 



is equivalent to a finite- dimensional symplectic map Hx • R 
{p',q'), which is called the constrained Toda (c-Toda) flow: 



R^\{p,q) 



p =9q, 

Aq~p-<q,q>q 

a 



q 



(5.11) 



Proof Directly calculate, and readily show ( ^.1C| ) <^=^ ( |5.11|) and {Hx)*{dpA 
dq) = dp A dq. I 
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When we understand the above two matrices {L'^'-')' and Mt in the following 



sense: (L^^)' L^^^, Mj 
variable), and set 



Tn U-e., q 



qn,p 



Pn, here n is the discrete 



Un = ±(< Aqn, Qn > - < Pn, qn > - < qn, qn >' 
Vn =< qn,qn >, 



(5.12) 



then the constrained Toda flow ( [5. 11] ) is none other than the well-known Toda 
spectral problem 

Llpn = {E~'^Un + Vn+ M„-E')V'„ = >^i^n, Efn = fn+1, E~'^ fn = fn-1 (5.13) 



with the above constraint ( p. 121) , A = Xj and ipn = qnj- Theorem 5.2 shows that the 
constrained Toda flow (Ht) has the discrete Lax representation (|5.10|) . Eq. ( ^.12|) 
is a kind of discrete Bargmann constraint of the Toda spectral problem ( [5.13| ). 

The Hamiltonian systems (Ht) and (Hq) share the common Lax matrix ( |5.1j ). 
Thus, they have the following same r-matrix: 



fx — X 



P-S, S 



f 1 \ 



-1 





(5.14) 



which is proven to satisfy the classical Yang-Baxter equation (YBE) 
[rij, rik] + [rij, rjk] + [vkj, r^fc] = 0, k = 1,2, 3. 



(5.15) 



5.2 The restricted Toda and CKdV flows 

Let us now consider the case on a symplectic manifold. We restrict the Toda and 
CKdV flows on the following symplectic submanifold M in i?^^ 



M = {{q,p) G R^^IF =<q,q> -1 = 0,G=< q,p> = 0}. 

Let us first introduce the Dirac bracket 

{/, 9}d = {f, 9} + ^({/, F}{G, g} - {/, G}{F, g}) 

which is easily proven to be a Poisson bracket on M. 

According to the thought of ref. [^, the following Lax matrix 



Ll"^ = Ll'^{X,p,q) 



i 



An{X) BniX) 
CniX) -AniX) 



(5.16) 



(5.17) 



(5.18) 
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yields 



- det Ll^ = Ix'TriLl^r = \x^+ < p, q > X + 2H^^ + ^ (5. 19) 
where 

= ^<Ap,q>-^<q,q><p,p>+^<p,q>^ (5.20) 

El^j = El^^{p,q)=p^q,-r„ i = l,...,N, (5.21) 

and Lq is defined by (5.2). 

An important observation is: if we consider the Hamiltonian canonical equa- 
tion produced by (|5.20|) in R^^, then this equation is exactly the well- 
known constrained AKNS flow, which will be discussed in the next subsection. 
Now, we first consider the Hamiltonian canonical equation restricted on M: 

{H^) : q, = {q, H^}d, p. = {p, H^}d, (5.22) 

which reads as the following finite dimensional system: 

Pa; = -|Ap+ i(< Aq,q > -l)p+ <p,p> q, 

q. = -p + lAq-l{<Aq,q>-l)q, (5.23) 
, < >= 1,< >= I- 

This is actually the CKdV spectral problem (^.7|) with the two constraints (Neumann- 
type) ig 

M=<Ag, g>-l, v=<p,p> (5.24) 

and A = Xj, ip = {pj,qj), j = 1,2, ...,N. So, the finite dimensional system ^5.23^) 
coincides with the restricted CKdV (r-CKdV) flow. 

Let us return to the Lax matrix ( [5.181 ). After endowing with an auxiliary matrix 
Mt,r as follows 

Mt,r=[\ (5-25) 

=< Aq — p, Aq — p > + < Aq, q > — < Aq, q >^, 

b =< Aq, q > -1, 

then we have the following theorem. 

Theorem 5.2.1 The discrete Lax equation 

(Ll'^yMT^R = Mt,rLI^, {Ll^y = Ll^{X,p', q') (5.26) 
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is equivalent to a discrete Neumann type of finite dimensional symplectic map Tir '■ 
{p, qf {p', q'Y 

p' = aq, 

q' = a-^{Aq-p-bq), (5.27) 
< <1,P >= h 
which is called the restricted Toda (r-Toda) flow. 

Remark 5.2.1 If we understand the above two matrices {L'^)' and Mt^r 
in the following sense: {^f)' {Ljf)n+i, Mt,r (MT,ij)„ (i.e., q ^ qn,P ^ 
Pn,ci an,b ^ bn, here n is the discrete variable), then the restricted Toda flow 
( p.27| ) on the symplectic submanifold M = {{q,p) G -R^^| < q,q >= l,< q,p >= |} 
is nothing but the discrete Neumann system studied by Ragnisco |]43[| . 

Let L^f = L^'"(A,p, g) ® / and L^f = I <^ L'j^{fi,p, q). Then, under the Dirac 
bracket ( [5.17| ) we obtain the following theorem. 

Theorem 5.2.2 The Lax matrix L'^^{X,p,q) defined by i\5.18i ) satisfies the 
following fundamental Dirac-Poisson bracket 

{Ll^{\) ? Ll^{^.)]n = h2(A,/.),LSf(A)] - hi(/i. A), L^^(/i)] (5.28) 

with a dynamical r -matrix 

ri2(A,/i) = P-5i2(A,/i), r2i(/i,A) = Pri2(^,A)P (5.29) 

ji — A 

where 



and P = |(/ + cTj (E) <7j) is the permutation matrix, dj (j = 1,2,3) are the 
Pauli matrices. 

This Theorem ensures that ( ^.21|) satisfies 

{EI%EI^^}d = 0, z,j,= 1,...,N. (5.31) 

For the r-Toda flow ([07| ), we have E^'^^ip', q') = Ej^^^ip, q) as well as YJ}=i Er^i =< 

p,q >= I from the discrete Lax equation ( p.26| ). Thus, in the set j-E^jj . ^, only 

E]^'^, -EJ*^, . . . , are independent on M. Therefore we obtain the following 

Theorem. 

Theorem 5.2.3 The restricted Toda flow Ht is completely integrable, and its 
independent and invariant (N-l)-involutive system is {E'^^}^Sl^ . 
For the restricted CKdV flow on M, we have 

^^ = ^EA,i^S (5.32) 
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which impUes 

{HgEl^^}D = 0, j = l,2,...,N. (5.33) 

Thus, the following Theorem holds. 

Theorem 5.2.4 The r-CKdV flow {H'^) is completely integrable, and its in- 
dependent (N-l)-involutive system is also {E]^f?^^Ii . 

Remark 5.2.3 As shown in this and last subsection, the r-Toda (i.e. Neumann- 
type) and the r-CKdV flows, and the c-Toda (i.e. Bargmann-type) and the c-CKdV 
flows respectively share the completely same Lax matrix, r-matrix and involutive 
conserved integrals. Thus, we say that the flnite dimensional integrable CKdV flow 
both restricted and constrained is the interpolating Hamiltonian flow of invariant of 
the corresponding Toda integrable symplectic map. 



5.3 The constrained AKNS and Dirac (D) flows 

Prom now on we assume: 



Lo = Lo(A,p,g) = f:-J— f ). (5.34) 



Let us again consider Eq. (5.18), and rewrite it as the following version: 



while we introduce 

= L^(A,p,?)=(^_°^ \^+L,. (5.36) 

Then we have 

1 N \ 2 pAKNS 

'-\'Tr{L^^^^)\\) =X' + 2X<p,q> + <p,q>' +2Hakns + E y ^ (5-37) 

1 1 ^ \^E^ 

-X^TriL^'fiX) = A2+A(< q,q> + <p,p>)--{<p,p> + <q,q >f+2HD+Y. 

(5.38) 



where 



Hakns = < ^p,q> ~ <q,q><p,p>, (5.39) 

Hd = ^{<Aq,q> + <Ap,p>) + ^{<p,q>^-<q,q><p,p>) 

+l{<p,p> + <q,q>?, (5.40) 

= 2pjq^ -T„ j = l,...,N, (5.41) 
Ef = p^^+q^-T^, j = l,...,N. (5.42) 
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Thus, Hakns and Hd generate the following two Hamiltonian systems 

m..s) : {'-^'^^r^^'''''^'" (5.43) 

i Px = =< P,P>q- Ap; 

^^^^ . I 9x = ^ =<P,g>(? + ^(<P,P> - < g,g>)p + Ap, ^^^^^ 
\ Px = = - <p,q> p-\[<P-,p> - <q,q>)q- ^q- 



It can be easily seen that {Hakns) and (i^i?) are changed to the well-known 
Zakharov-Shabat-AKNS spectral problem 



Vx 



A u 
V -A 



(5.45) 



and the Dirac spectral problem |30 



Vx 



X — u 



-X — u 



(5.46) 



with the constraints u = — < q,q >, v =< p,p >, A = A-,-, y = {qj,pj)^, and the 
constraints u = — 1(< p,p > — < q,q >), v = — < p, q >, X = Xj, y = {qj,Pj)'^, 
respectively. 

Therefore (Hakns) and {H^,) coincide with the constrained AKNS (c-AKNS) 
system and the constrained Dirac (c-D) system, respectively. 

Let Lf{X) = L^(A) ® /, L^(/x) = / ® L^(/i) (J = AKNS, D). Then we have the 
following theorem. 

Theorem 5.3.1 The Lax matrices L-^{X) (J = AKNS,D) defined by Eq. 
( ]5. 33^ ) and Eq. ( |5. satisfy the fundamental Poisson bracket 



{L-'iX) ? L\^^)} = [n2(A,/i),Lf(A)] - hi(/x,A),L2^(/i)]. 



(5.47) 



Here the r-matrices ri2(A,/i), r2i(/i. A) are exactly given by the following standard 
r -matrix 



n2(A,/i) 



/i — A 



P, r2i{fi,X) = Pr,2{fi,X)P, 



(5.4J 



P 



/ 1 \ 

10 

10 

1 



1 ^ 

^ i=l 



CTi). 



(5.49) 



So, the c-AKNS and c-D flows share the same standard r-matrix ( ^.48|) , which 
is obviously non-dynamical. However, the two constrained flows, produced by 
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(|5.45|) 's and ( |5.46| )'s extensive spectral problems ( |6.1|) and ( |6.2|) (they are 



gauge equivalent), have different r-matrices (see section 6). 

Remark 5.3.1 In fact, the r-matrix ri2(A,/i) in the case of the c-AKNS and 
c-D flows can be also chosen as 

ri2(A,/i) = ^-P + /®^, 5= f M (5.50) 
yU — A \ c a J 

where the elements a, b, c, d can be arbitrary C°°-functions a(A, ^,p, q), b{X, fi,p, q), 
c(A, fi,p, q), d{X, fj,,p, q) with respect to the spectral parametres A, /i and the dynam- 
ical variables p, q. This shows that for a given Lax matrix, the associated r-matrix 
is not uniquely defined (there are even infinitely many r-matrices possible). Here 
we give the simplest case: a = 6 = c = c? = 0, i. e. the standard r-matrix (|5.48|) . 

5.4 The constrained Harry-Dym (HD) and Heisenberg spin 
chain (HSC) flows 

The constrained Harry-Dym system describes the geodesic flow on an ellipsoid and 
shares the same r-matrix with the constrained Heisenberg spin chain (HSC). To 
prove this, we consider the following Lax matrices: 



= L--(A,p,,)=f A-+<^>^>A- 1+^^^ (5.51) 

< p,p > X < p,q > X ' 



\ — < p,p > \ < p,q > A J 

Here L^^^ is included in the generalized Lax matrix (2.1), but L^^ is not. We 
need two associated auxiliary matrices 

Mhd = f <Ap,p> , ), (5.53) 

(-iX<Ap,q> iX<Aq,q>\ , . 

Mhsc = -x a -x a (5-54) 
\ —iX < Ap,p > tX < Ap, q > J 

Theorem 5.4.1 The Lax representations 

L^"" = [Mhd,L''% (5.55) 

L^'^ = [Mhsc,L^'^) (5.56) 
respectively give the following fininte dimensional Hamiltonian flows: 

q.=P=^\TQ--^, 

(Hhd) : { P. = -ii^Aq = -^Itq-^, (5.57) 
< Aq,q>= 1; 
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■^^^ ■ = i < Ap,q > Ap - i < Ap,p > Aq = ' 
wt/i t/ie Hamiltonian functions 

i^/fD = -<p,p>--—^^^{<Aq,q>-l), (5.59) 
2 2 < A^g, q > 

Hhsc = < Ap,p>< Aq,q> -^i < Ap,q>^ . (5.60) 
In Eg. ( |5.5?| j TQ^~^ is a tangent bundle in R^^ : 

TQ^-^ = {{p, q) e R^^\F =< Aq, q> -1 = 0,G =< Ap, q >= 0}. (5.61) 
Obviously, Eq. ( |5.57| ) is equivalent to 

q.. + ^4|l^Ag = 0, < Aq, q >= 1, (5.62) 
< A^q, q > 



which is nothing but the equation of the geodesic flow on the surface < Ag, q > = 
1 in the space and also coincides with the constrained HD (c-HD) flow [^. In 
addition, Eq. ( p.58| ) becomes the Heisenberg spin chain spectral problem 



-iXw —iXu 
-iXv iXw 



y, t' = -l, (5.63) 



with the constraints u = — < Aq,q >, v =< Ap,p >, w = — < Ap,q >, A = 
■^i? y = {qjiPjY ■ Thus, Eq. ( [5.58[ ) reads as the constrained Heisenberg spin chain 

(c-Hsc) flow gg. 

Their Lax matrices (|5.51|) and ( |5.52|) share all elements except one, namely 

(:':)■ 

This element does not affect the calculations concerning the fundamental Poisson 
bracket, one can readily deduce that the c-HD flow and the c-HSC flow possess the 
same non-dynamical r-matrix 

2A 

ri2(A, li) = ^^^_^^ P, r2i(/x. A) = Pri2(/i, A)P. (5.64) 

Remark 5.4.1 The r-matrix ( ^.64|) of the c-HD and c-HSC flows can be also 
chosen as 

2A 

ri2(A, ^l) = -. -^P + 10 S. (5.65) 

fi{fi - A) 

Evidently, Eq. ( p.64| ) is the simplest case: S* = of Eq. ( p.65| ). 
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5.5 The constrained G and Q flows 

In this subsection, we introduce the following Lax matrices: 

L = L (A,p,,)=(^ ^ -(l+<p,,»A- J+^°' ^'-''^ 

L« = L«(A,p,,) = (-^^'' <''^i;^"')+i^o. (5.67) 
If we set 

/ -(< p, n > - < (7, (7 >) - 1 \ 

Mg = 1 , " " 1 , , 5.68 

/A + ^<A,,g><p,p> ^<A.,.> I 

\ -^<p,p>A -A - 2^ < Ag,g ><p,p > ' 

with 



a = ^J{< p,p > - < Aq,q >y - A < Aq,p >, l3 = l-<p,q>, (5.70) 

then, by a lengthy and straightforward calculation we obtain the following theorem. 
Theorem 5.5.1 The following Lax representations 

L^=[Ma.L''] (5.71) 

and 

L^ = [Mq,L'^] (5.72) 

where the first one is restricted to the surface Mi = {{p,q) G -R^^| < p,q >= 
0, < Aq,q >< p,p > + < Aq,p >= 0} in the space B?^ , respectively produce the 
finite- dimensional systems: 



and 



qx = i(-Ag + {<p,p> - < Aq,q >)p) -p, 
p^ = MAp+{<p,p > - < Aq,q>)Aq)+Aq, 



q^ = Aq+^ < Aq,q> p+^ <p,p>< Aq,q> q, 
Pec = -Ap - ^ < p,p > Aq - ^ < p,p >< Aq,q > p. 



(5.73) 



(5.74) 



Eqs. (|5.73|) and ( |5.74| ) turn out to be the spectral problem studied by Geng 



(simply called G-spectral problem 

X{v+1) Xu 



y^-( .r"\ '":']y (5.75) 
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with the constraint condition 
1 



1 



u 



^ yi<p,p>-<Aq,q>y-4:<Aq,p> 

< p,p > — < Ag, q > < p,p > — < Aq, q > 



a 



< p,p > — < Aq, q >y — 4 < Aq,p > 



A = \j, y = {qj,Pj)'^, and the spectral problem proposed by Qiao (simply called 



Q-spectral problem) 



A — ^uv 
Xv 



u 

-A + \uv 



(5.76) 



with the constraint condition 



u 



< Aq, q > < Aq, q > 



< p,p > 



1— < q,p > 
< p,p > 

1— < q,p > 



X = Xj, y = {qj,pjY , respectively. 

So, Eqs. ( |5.73|) and ( ^.74| ) are the constrained Geng (c-G) flow and the con- 
strained Qiao (c-Q) flow, and they have the same non-dynamical r-matrix: 



n2(A,/i) 



2 2 

-P - -S, S 

IX — X /i 



/ \ 



10 





(J_ (X) (J 



(5.77) 



Here, the r-matrix ri2(A,/i) can be also chosen as 

'^12(A,/i) 



-^P --S + I ®S. 
fi — X fi 



(5.78) 



Eq. (15.771) is the simplest case: 5* = of Eq. (|5.78|) . 



We have already seen that the r-matrix ri2(A,/i) satisfying the fundamental 
Poisson bracket is not unique (in fact, infinitely many) and is usually composed of 
two parts, the first one being their main term, and the second one being the common 
term I ^ S. Usually, to prove the integrability we choose their main term as the 
simplest r-matrix. 
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6 An equivalent pair with different r- mat rices 

This section reveals the following interesting fact: a pair of constrained systems, 
produced by two gauge equivalent spectral problems, possesses different r-matrices. 
In 1992, Geng introduced the following spectral problem EH 



M0, M 



iX — ijSuv u 

V —iX + i[3uv 



-1 



(6.1) 



where u and v are two scalar potentials, A is a spectral parameter and /5 is a con- 
stant, and discussed its evolution equations and Hamiltonian structure. Eq. ( |6.1| ) is 
apparently an extension of the AKNS spectral problem ( |5.45| ). Two years later the 
author considered an extension of the Dirac spectral problem (|5.46| ) 



— IS 



A + r + /5(s2 - 
is 



(6.2) 



^ -A + r - /3(s2 - r^) 
where r, s are two potentials, and obtained a finite dimensional involutive system 



being not equivalent to that one in ref. |2^. But, the spectral problems (^TT]) and 
are gauge equivalent via the following transformation 



i) = G(t), G 



1 1 

i —i 



(6.3) 



V = i{r — s), u = —i{r + s). In ref. |Q, Wadati and Sogo discussed the gauge 
transformations of some spectral problems like Eq. |5.63| . 



Now, we discuss their r-matrices. Let us consider the following two Lax matrices 



rGX 



rGX 



iQZ ^ LQZ(^x^p^q) 







1 + 2i/3 < p, g > 






-1 - 2il3 < p,q > 



■ /3(< p,p> + < q,q>) 



+ f3{<p,p> + < g,g >) 

Then calculating their determinants leads to the following Hamiltonian systems 

<p,p><q,q> ^ <Q,q> 



(6.4) 
(6.5) 



dp 



■2i/3 



aq r (l+2if)<p,q>)'^ l+2ip ^ ' 



(6.6) 



and 



qx 

Px 



dp 

OHqz 

dq 



An _ /3 4<p,g>^+(<p,p>-<g,g>)^ _ 
^^P P (l-2/3(<p,p>+<g,g>))2 P 

^2 I „^„^„ „■^^2 



2<p,g>g+(<p,p> — <g,g>)p 
l-2/3(<g,g>+<p,p>) ' 
A „ 1 o 4<p,g>2 + (<p,p>-<g,g>)2 , 2<p,q>p-{<p,p>-<q,q>)q 
IVq^p (l_2/3(<p,p>+<g,g>))2 q^ i_ 



-2/3(<g,g>+<p,p>) 



with the Hamiltonian functions 

Hgx = i< Aq,p> 



< p,p >< q,q > 
' 2(1 + 2ip <p,q>] 



(6.7) 



(6.8) 
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and 



QZ = - < Ap,p> +- < Aq,q> —- — — — . 6.9 

2 2 4 — 8p[< p,p > + < q,q >) 



Obviously, Eqs. (|6^ ) and (|6.71 ) become Eqs. (|6.1|) and 

< q,q > < p,p > 



u 



1 + 2ii3 < p,q >' 



1 + 2iP < p,q > 



with the constrants 
(6.10) 



A = \j, = {qj,Pj)^, j = I, . . . , N; and the constraints 



-2i < p,q > 



- < p,p > + < q,q > 



(6.11) 



1 - 2f3{< q,q> + <p,p>)' 1 - 2/3(< q,q> + <p,p>)' 
X = \j, ip = {qj,Pj)^, j = l,...,N, respectively. Thus, the finite dimensional 



Hamiltonian systems ( |6.6| ) and ( |6.7| ) are respectively the constrained fiows of the 
spectral problems ( |6n|) and ( |6.2|) . Since they have Lax matrices (|6.4| ) and (|6.5| ), 
then the r-matrices of (|6.6|) and ( |6.7|) are respectively: 



ri2(A,/i) 



/i — A 



P + AipS, S 



and 



n2(A,/i) 



/i — A 

which are apparently different. 



P + 2pS, S 





( 1 














































lo 












( 











-1 


\ 










1 












1 










\ 


-1 











/ 



(6.12) 



(6.13) 



7 New integrable systems 



In this section, three new integrable systems are generated as the representatives 
from our generalized r-matrix structure. 

1. The first system is given by case 6 in section 4. The corresponding r-matrix 
and involutive systems are respectively 



'^12(A,/i) 



fi — X 



P + -S, s 



/ 1 \ 



10 

1 



(7.1) 
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and 



Ej = 2(< p,q> +c)Xj ^pjqj+ <q,q> ^p] - T^, j = 1, . . . , AT, 



(7.2) 



where c & R. Thus, the finite dimesional Hamiltonian systems {F^ defined by 
i^i. = Ef=iAfi?j,m = 0,...,i.e. 

Fl = 2{<p,q>+c)<X'p-'p,q> + <q,q><X'P''p,p> 

- {< A'q,q >< A.^p,p > - < A'q,p >< A.^p,q >) (7.3) 

i+j=m—l 

are completely integrable. Particularly, with m = 2 the Hamiltonian system (-^2^): 

dFi 



Qx = ^ = 2cAg -2 < Aq,q> p + A< Ap,q> q + 4:<p,q> Aq, 
Px = ~~Q^ = —2cAp + 2<p,p>Aq — A< Ap, q > p — A < p,q > Ap, 

is a new integrable system, which becomes the following spectral problem 

('(2c + 4«)A + 4,. -2.. \ 

^ I 2s\ -{2c + Av)X-Au ; ^ 



(7.4) 



(7.5) 



with the constraint conditions u =< Ap,q >, v =< p,q >, w =< Aq,q >, s =< 
p,p >, and A = \j, = {qj,pj)^, j = 1,. . . ,N. Apparently, the spectral problem 
( [7.5|) is new. 

2. The second system is produced by case (7.1) in section 4. The correspond- 
ing r-matrix and involutive systems are respectively 

/ g'2 o\ 

9 9 f 

n2(A,/x) 



2 2 

-P + -S, S 

ji — X n 



f'2 



-1 



V 



0-100 




(7.6) 



and 



= 2c\fp,qj + {<q,q> +g2)\fp) + {<P,P> -f2)Xj\] - T,, j = 1, . . . , iV (7.7) 

where c e -R, /2 = /2(< P,q >), 92 = 92{< Vil >) ^ C'^{R). Hence, the 
Hamiltonian system (^2^) defined by F^ = J2jLi ^'jEj, i.e. 

= 2c<Ap,q> +2< Ap,q><p,q> +gi<p,q>) < Ap,p > -f{<p,q>) <Aq,q> 

(7.8) 

is completely integrable. Meanwhile the Hamiltonian system (-Ff): 



Px 



.211 

dq 



2cAq + 2 < Aq,p> q + 2 <p,q> Aq 

+2g{< p, q >)Ap+ <Ap,p> g'{< p, q >)q- <Aq,q> /'(< p, q >)q, (7.9) 
-2cAp - 2 < Aq, p>p-2<p,q>Ap 

- <Ap,p> g'{<p,q>)p + 2f{<p,q>)Aq- < Aq, q > f {< p, q >)p, (7.10) 
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can also related to a new 2x2 spectral problem with some constraint conditions. 

3. The third system is derived by case (7.4) in section 4. The corresponding 
r-matrix and involutive systems are respectively 



n2(A,/i) 



2 2 

II — A 



/ 1 \ 



10 

\0 c'_^ 1 J 



(7.11) 



and 



(7.12) 



E] = 2{- <p,q> +c)Xj ^Pjqj+ <q,q> \ ^v] " C-iA^- - 

where c E R and c_i = c_i(< p, g >) G C'^{R). The following Hamiltonian system 
{Fi): 



1^ = ^ = 2cA9- <Ag,q> (c'_i(< p, q >)q + 2p), 



Pa; 



dp 
OF: 



= -2cAp + 2(c_i+ < p,p >)Aq + c_i(< p,g >) < Ag, g > p, 



(7.13) 



is one of their products, where 

F2 = 2c < Ap,q > - < Aq,q > {c^i+ < p,p >). 



(7.14) 



In general, with any c_i Eq. ( [7.13|) can't be changed to a 2 x 2 spectral problem 
with some constraints. But with two special c_i: c_i = and c_i =< p,q >, Eq. 
( [7.13| ) can respectively become the spectral problem |^ 



2cA -2v 
2u\ -2cA 



(7.15) 



with the constraint conditions u =< p,p >, v =< Aq, q >, and the spectral problem 

(7.16) 



2cA - V -2v 
2uX -2cA + V 



with the constraint conditions u =< p,q + p >, v =< Aq, q >. Here in Eqs. ( |7.15 ) 
and ( |7.16| ) A = Aj, (j) = {qj,pj)'^, j = 1, . . . , N are set. Eq. ( |7.16| ) is a new spectral 
problem. 

Remark 7.1 We can consider further new integrable systems induced by 
Theorem 3.1. 

Remark 7.2 The above procedure actually give an approach how to connect 
an r-matrix of finite dimensional system with a spectral problem, which is closely 
associated with integrable NLEEs. 
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8 Algebro-geometric solutions 



The ideal aim for nonlinear differential equations is of course to obtain their explicit 
solution. In this section, we connect the finite dimesional integrable flows with inte- 
grable NLEEs, and solve them with an explicit form of algebro-geometric solutions. 
Here, we take two examples: one being the periodic or infinite Toda lattice equa- 
tion, the other the AKNS equation with the condition of decay at infinity or periodic 
boundary. 



8.1 Toda lattice equation 

The Toda hierarchy associated with Eq. ( p.l3| ) is derived as follows: 



V.n. 



JG"^, J = 0,1,2,... (8.1) 



where {G" = J ^KG"_i}°^q is the Lenard sequence, G^^ = G KerJ, for 

all a = a(tj), [3 = Pitj) G C°°{R), the two symmetric operators K, J are 



^ ^ ( \un{E - E-^)Un Un{E-l)Vn\ Un{E - 1) 

\vn{l-E~^)Un 2{ulE - E-^ul) ) ' \{l-E~^)Un 



i2) 



In particular, with j = 0, /3 = 1 Eq. reads as the Toda lattice 

Un = Un{Vn+l-Vn), ?)„ = 2(m^ - J (8.3) 

which can be changed to 

x„ = 2(e2(^"+i-^") - e2(^"-^"-i)) (8.4) 
via the following transformation 

Un = 6^"+^-'=", Vn = Xn- (8.5) 

It is easy to prove the following theorem. 
Theorem 8.1.1 

1) Let Gn = {G^,^\G^^^f, VGW,^^^) G C^{R). Then the operator equation 
ViGn), L] = L,{KGn) - L,{JGn)L 
possesses the operator solution 

V{Gn) = -{E-\^)G^^^E-^ + ^((E-i«„G«) - «.G«) + u^G^^^E (8.6) 
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where [■, ■] is the usual commutator; the operator L is defined by Eq. ( ^.ISj ); = 
E-%+^2 + ^iE, Ve = (^1,6)^,6,6 e C°°{R). 

2) Let us choose the special Gn = G^, j = —1,0, 1, then the Toda hierarchy 
( \8.1{ ) has the following Lax representation of operator form 

Lt^ = [WiG]),L], J =0,1,2... (8.7) 

where the operator W{G]) = Ei=o '^(G'Li)^^'"''- 

Particularly, the standard Toda equation ( p.4|) possesses the Lax representa- 
tion of operator form Lt = [W{Gq), L], where the operator W{Gq) = e^"+^~^"E — 
^x„-x„-i ^ and m„, Vn in L are substituted by 

We have shown that the c-Toda flow and the c-CKdV flow share a common 
nondynamical r-matrix and in particular, this ensures the integrability of their flows. 
A calculation of determinant yields their common A^-involutive systems 

E^ = XM-pl-<P,q>ql- E a = 1, . . . , iV (8.8) 

which are independent and invariant (i.e. Ea{\,p,q) = Ea{\,p' ,q'))- Apparently, 
the functions Eg = J2a=i ^a^a, s = 0, 1, 2..., are given by 

Fs = < A''+^p,q > - < A''p,p > - < p,q >< A''q,q > 

- {< A^p,p >< A''q,q > - < A^p,q X A''q,p >) (8.9) 

j+fe=s— 1 

and {Ejn, Ei} = 0, Vm, / G which implies the Hamiltonian systems (Eg) are 
completely integrable. 

Let {po(ts),qo(ts))^ be a solution of the initial problem 



— ( ^ \ = ( \ I P\ =( 

dts[qj [dEjdp J' UA^=o Uo 



(8.10) 



Set 



f 'ii )=H}( ^fi ) (8.11) 

where Ht is defined by Eq. ( |5.11| ). Now, we rewrite Eq. ( [5. 121 ) as a map / : R^^ — *■ 
i?^ defined by 

/: (Pn,gn)'^l >{Un,VnV. (8.12) 

Then, we have the following theorem. 

Theorem 8.1.2 {un(ts),Vnits))'^ = f{pn(ts),qn(ts)) satisfies the Toda hierar- 
chy 

d ( Ur,.^ 



dt. \ v„ 



JG'l, s = 0,l,.... (8.13) 
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Particularly, with s = the following calculable method 

\qo J V (io{t) J V Quit) J V Mt) J 



produces a solution of Toda lattice equation (|87^ ). Thus, the standard Toda equation 
j4[) has the following formal solution 

Xnit) = f < qn{t),qnit) > dt. (8.15) 



We shall concretely give the expression < qnit),qn(t) >. 
Let us rewrite the element CtcW of Eq. ( p.l| ) as 

^^^^^^^-^V K{X)=m^-K), (8.16) 
and choose distinct real zero points /ii, . . . , fijy of Q{X). Then, we have 

N N N 

Qi\) = l[i\-fi,), <q,q>=J2K-J2fx,. (8.17) 

j=l a=l j=l 

Let 

TT, = Arcif^j), (8.18) 

then it is easy to prove the following proposition. 
Proposition 8.1.1 

{Hi,fij} = {7ri,nj} = 0, {TCj^fii} = 5ij, j = 1, 2, A^, (8.19) 

i.e. TTjyfJ'j are conjugated, and thus they are the seperated variables /^^. 

Write detL^^(A) = -A^^W - Btc{\)Ctc{\) = - ELi = "TO' 

where is defined by (8.8), and -P(A) is an + 2 order polynomial of A whose first 
term's coefficient is |, then 7r| = j = 1, ■ ■ ■ ,N. Now, we choose the generating 

function 



N N 



W = t W,{f,„ {E^}^^,) = t r^"^ g^dX (8.20) 
j=i j=i-^Hio) \K{X) 

where /ij(0) is an arbitrary given constant. Let us view Ea (a = 1, . . . , A^) as actional 
variables, then angle-coordinates Qa are chosen as 

dW 

Qa = wTr, a = l,...,N 

OJ^n 
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i.e. 



<5q = V / Wa, = — / ^dA, a = l,...,A^. (8.21) 

fc=i-^/^fc(o) 2Ji^(A)P(A) 



Hence, on the symplectic manifold (i?^^, c/E'a A (iQa) the Hamiltonian function 
Fq = J2a=i Fa produces a hnearized flow 



as,' 



Ea = 0, 



thus 



f^J.k{n+l) ~, 

Ejn)=EJn-l] 



Qa (n) = go + t + C„n, C« = ELi Qn) ' ' '^'^ 



(8.22) 



(8.23) 



where Ca are dependent on actional variables {Ea}a=i, and independent of t; (5° is 
an arbitrary fixed constant. 

Choose a basic system of closed paths ai, Pi, i = 1, . . . , N of Riemann surface 
T: fj? = P{X)K{X) with N handles, ujj {j = 1,...,N) are exactly N linearly 
independent holomorphic differentials of the first kind on this Riemann surface T. 
LVj are normalized as tOj = J^iLi^j^ii^i, i.e. tOj satisfy 

tOj = Sij, <j>^ U!j = Bij 

where B = {Bij)]^^.^ is symmetric and the imaginary part ImB of -B is a positive 
definite matrix. 

By Riemann Theorem ||25| we know: /ifc(n) satisfies J2k=i I^^{o) ~ ^i' ~ 
(f)j{n, t) = Y^iLi '^iAQ^ + t + cin), j = 1, . . . , iff /ifc(n) are the zero points of the 
Riemann-Theta function 0(P) = Q{A{P) — cf) — K) which has exactly zero points, 
where A{P) = {jj^^uj,,- ■ ■ , /^^ umY, = <P{n, t) = (0i(n, t), ■ ■ ■ , 0^(n, t))^, G 
is the Riemann constant vector, Pq is an arbitrary given point on Riemann surface 
f. 



Because of 18 



— / Arflne(P) = Ci(r) (8.24) 

2TXi J-y 

where the constant Ci(r) has nothing to do with 0; 7 is the boundary of simple 
connected domain obtained through cutting the Riemann surface T along closed 
paths ai. Pi. Thus, we have a key equality 

N _ 

J2fik{n) = Ci{T) - Resx=oc,Xd\ne{P) - Resx=oc,Xd\ne{P) (8.25) 

k=l 
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where ooi := (0, ^P(^-i)i^(^-i)|,=o), 002 := (0, -^P(z-i)ir(^-i)U=o). Through 
a lengthy careful calculation and combining (|8.17|) , we obtain 

, ^ d QUinA) + K + rii)\ 

< >= E A. - C,{T) + - In (8-26) 



dt\ Q{(f){n,t) + K + 7]2) 



where the j-th component of rji {i = 1,2) is r/j.-,- = J^-UJj. By the Riemann surface 
properties, we can also have Y.fLi '^j,ici = j^l ujj = J2iLi "^3,1 which imphes 



Q = J^^ Coi = Jp^' l^^^^=^dX. So, the standard Toda equation (|8j) has the 
following explicit solution, called algebro-geometric solution 

= e(Xt)?F/?z) ^ ^ « ^ ''-''^ 

where U = RC, V = RJ, Z = RQ"^ + K + 7]i with C = (ci, . . . ,C7v)^, J = 
(1, . . . , 1)^, QO = (QO, . . . , Q^)^, P = Ell K - C,{T), while matrix R = (r,, ^^xiv 
is determined by the relation YA=i^j,i§ai^i — ^iji ^^"^ ^ certain constant which 
can be determined by the algebro-geometric properties on the Riemann surface 
f The symmetric matrix B = (-Btj)ivxiv in function is determined by 



I]- 



Hence, the algebro-geometric solution of Toda lattice (^]3|) is 

^ „x„+i-x'„ ^ . e^{U{n+l)+Vt+Z) 

^ ^ e{U(n+2)+Vt+Z)e{Un+Vt+Z)'' (0.00^ 

-Pi In ^jUri+Vt+Z) l^'^^J 
■^n dt 0(f/(n+l)+Vt+Z) ■ 

Obviously, the algebro-geometric solution Unit) and f„(t) given by ( ^.28| ) are quasi- 
periodic functions, and they are periodic iff f/ = where M is a A^-dimensional 




integer column vector. It is easy to see that (^.28) is the finite-band solution of Toda 



lattice (|8.3| ) if Ai, . . . , A^v are chosen as the eigenvalues of Toda spectral problem 

(EH)- 



8.2 AKNS equation 

In subsection 5.3 we have shown that the constrained AKNS flow shares a common 
r-matrix with the constrained Dirac flow, therefore they are integrable. Now, we 
derive the algebro-geometric solution for the second order AKNS equation (|8.31|) . 
It is well-known that the AKNS hierarchy is given by 

= JG„ i = 0,1,2,... (8.29) 
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where {Gj = J ^KGj-i}'^^^ is the Lenard sequence, G_i = (0,0)-'" G KerJ, the 
two symmetric operators K, J are 



I 2.a-„ ^-^''fM, J = 2f° -1 1. (8.30) 



A representative equation (j = 2) of ( ^.291) is 



Ut = -^U^x + U^V, Vt = ^Vxx - v^u, t = t2. (8.31) 

The independent A^-involutive system of the constrained AKNS flow is expressed 
by Eq. ( |5.41| ). Similarly, we consider the following Hamiltonian functions 



N 



3 = ^ 



= 2<A''p, g>- ^ {< A^p,p >< A''q,q > ~ < A^p,q >< A''q,p >). (8.32) 

j+fc=s-l 

Let {p{x,ts),q{x,ts))'^ be the involutive solution of the consistent Hamiltonian 
canonical equations (Hakns), {F^^'^^). Then, we have the following theorem. 

Theorem 8.2.1 u = — < q{x,tj),q{x,tj) >, v =< p{x,tj),p{x,tj) >, j = 
0,1,2..., satisfy the higher-order AKNS equations ^8.2^^ ). Particularly, Eq. ^S.St^ ) 
is solved with the following solution: 

u = - < q{x,t2),q{x,t2) >, V =< p{x,t2),p{x,t2) >, (8.33) 

where {p{x,t2),q{x,t2))^ is the involutive solution of the consistent Hamiltonian 
systems {Hakns), (^2^^^^). 

In the following procedure we shall express Eq. ( |8.33| ) in an explicit form of 
algebro-geometric solution. To do so, let us rewrite Eq. (|5.35|) as follows: 

j^AKNS ^ / AakNs{.>^) BaKNsW \ g^-j 



where 



N ^ 

AaknsW = 1 + rPj^J^ (8-35) 

N ^ 

BaknsW = — r^i' (8-36) 
j=l ^ - 

N ^ 

CaknsW = El-^P?- (8-37) 

j=i ^ - 
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BaknsW, CaknsW can be changed to the following fractional 



form: 



R ... <q,q>QBiX) ^ ... _ 
^AKNSK'^) = T77V\ ' ^AXAfsl^J = 



< 



p,p > QcW 

K{X) 



(8.38] 



where 



N N 

<q,q>QB{\) = n (A -A,), 

j=l k=l,k^j 
N N 

<p,p>Qc{\) = Y.P] n (A -A,), 

o — 1 i. — 1 h-A^ 



K{\) 



j = l k = l,kytj 

N 

n(A-A,). 



Respectively choosing — 1 distinct real zero points /if, ... , fi^-i and fii 
of QbW and QcW leads to 



ti ti <p,p> 

N-i N-1 < n > 

1,-1 1,-1 ^ q^ q ^ 



k=l 
N-l 



k=l 



, . . . , fJ^N-l 



8.39) 
8.40) 
8.41) 



N-i N-l < A^n n > 

^1 t^i <p,p> 

where Ai = Y.f=i Aj, A2 = Y.k,j=i,j<k AjAfc are two constants. One hand, Ux = —2 < 
q, q, >= -2 < q, ^^^^ >=' -2 < Ag, g > -2t 
fixed function of t. Thus from Eq. ( p.40| ) we have 

d 



2 < Ag, g > —2uco(t), here Co(t) is an arbitrarily 
l- ( p.40| ) we have 

— lnu = 2Ai-2E/^f-2co(t). (8.42) 

ox f-; 

>= -2 < A2g,g >. This 



, ut, = -2 < g,gt, >= -2 < g, >= - ^2 



On the other hand, - - 1711-2 - - - 

is combined with Eq. ( ^.40|) to give the equality 



o N-l N-l 

^ ln« = (A, - E f^^? - E (-"f - 2^2 + 

•^^2 k=l 



Af-1 

E 

fc=i 



(8.43) 



So, we obtain 



M(x,t) = u{xo,to)exp{ [ 

J to 



N-l 



+ 



iV — i iV — i 

(^1 - E /^f - E (/^f - 2^ + ^? 

fc=l k=l 
N-l 

2Ai - 2 E /if - 2co(t) 
fc=i 



dt 



dx), t = t2, (8.44) 
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where xo, to are two fixed initial values. Similarly, v{x,t) has the following repre- 
sentation 

iV-l N-1 

{Al - E Z^?)^ - E i^^k? - 2A, + Al 

*o L k=l k=l 

N-1 



v{x,t) = f (xo, to) exp(- 



dt 



2Ai - 2 ^ /i^ - 2co(t) 



k=l 



dx), t = t2 



(8.45) 



Since Eqs. ( p.44| ) and ( ^.45| ) solves nonlinear equation ( p.31| ), then in order to 
obtain their explicit form it needs calculating the four key expressions Y^k=i{l^kYi 
J = B,C] k = 1,2. For that purpose, we follow the approach in the case of Toda 
lattice equation. For the present two set of Darboux coordinates /i/, J = B,C; 
J = 1, 



— 1, we have the following key equalities like Eq. (^.251) 

N-1 2 

^ (/i/)'= = CkiT) - J2 Resx=oo^'d\n e(A(P) -<P-Kj), 

j=l s=l 

J = B,C; k = l,...,iV- 1, 



(8.46) 



where Ck{T) is a constant [p6| , |56| only determined by the compact Riemann surface 



r {genus = N 
(A - Afc); ooi 
A{P) 



-1): fi ^ = Pakns{\)K{\) , PaknsW = i^(A)+E.f=i 

: (0, VPaKNs{z-')K{z'^)\,=o), CX)2 



r^AKNS 



n 



TV 

kj^j,k=l 



(0, -VPAKNs{z-')K{z-^)U=oy, 



Jp^^ 00 is an Abel map in which Pq is an arbitrarily fixed point on F, uj 



[UJi, UJN-1) 



El 



N-1 



n 



JV 



.i(A-Afc) 



dX is a normalized holo- 



'-^ ^'^2y/K{\)PAKNsW^ 

morphic differential form, and rj^i is the normalized factor; The j-th component 
(pj{x,t) oi N — 1 dimensional vector equals to J2hJi^ ''"jAQ^ + ^^i^ + ^■^f'^ + 
Ci(t) + Ci{x)) with the arbitrary constant (5° and functions Ci(t), Ci{x) G C°°{R); 
Kb, Kc G C^~^ are the two Riemann constant vectors respectively associated with 
the Darboux coordinates /i^, /i^; Riemann- Theta function ||3^ ©(0 is defined on 
Riemann surface F. 



A lengthy calculation of Residue at oo^, s 



N-1 



Ci(r) 



N-1 

E(/^, 



— In 

dx 69 ' 



= 1,2 for k 



1, 2 yields 



J\2 

j > 



dx'^ 



In 85^0^ 



(8.47) 
(8.48) 



where 6f 
of the N - 



= e{(l) + Kj + ris), J = B,C, 7]s,j 
1 dimensional vector rjs- 



2 



1, 2 is the j-th component 



Substituting the above equalities into (|8.44|) and ( ^.45| ), and sorting them, we 
obtain the explicit solution of nonlinear equation ( |8.31| ): 

(x-xo)^l 



u{x, t) 



u{xo,to)d 



a(t-to)+2{b-co(t)){i 



^2 \2\ 



QB\t=toiQB) 



\X=XQ 
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ef 

x^exp, 



to 



f)2 f) 

— In e^G^ + ffc + — In -^)^ 



dt), 



.49) 



v{x, t) 



v(rn f^)p-<^it-to)-2{b-co{mx-xo)'^\ 



X 



^2 



c exp( 



to 



dt), (8.50) 



where a = Aj - C2{T) - 2^2, b = Ai - Ci{T) are two constants, co(t) G C°°(i?) is 
an arbitrarily given function of t, and xq, to are the initial values. Therefore, we 
have the following theorem. 

Theorem 8.2.2 The AKNS equation (^. 31\ ) has the explicit solution l \8.49l ) 
and ^8. 5C\ ) given by the form of Riemann-Theta function, which is called the algebro- 
geometric solution. 

An analogous calculational process will lead to the algebro-geometric solution of 
the higher-order AKNS equation (|8.29|) . 



9 Conclusions and problems 

In finite dimensional case. Lax matrix is enough to provide many important in- 
tegrable properties like r-marix, Hamiltonian, integrability, Darboux coordinates, 
and even later algebro-geometric solution. Therefore we specially stress to use Lax 
matrix instead of Lax pair in finite dimensional case. 

The generalized r-matrix structure is given to emphasize the classification and 
united sketch of finite dimensional integrable systems. We have already seen that 
only is there one concrete r-matrix structure, then the corresponding Hamiltonian 
flows are surely integrable and even in some cases the associated spectral problems 
are new. 

In the paper, we develope our generalized structure to solve some integrable 
equations with algebro-geometric solution. This is an extension of nonlinearization 
methods It is found that this procedure can be also applied into other integrable 
NLEEs [p6| p5| , p^ . In this sense, we successfully realize a procedure from finite 



dimensional flows to infinite dimensional systems when we have some constrained or 
restricted relation between them. Of course, there are still other methods to solve 



integrable NLEEs. Recently, Deift, Its and Zhou [0, |T^ obtained the 0-function 
solutions of some integrable NLEEs like the KdV, MKdV, nonlinear Schrodinger 
equation by using Riemann-Hilbert asymptotic method. All these methods are still 
under the development. 

It should be pointed out that our procedure is carried in the symplectic space 
{R'^^ , dp A dq) (i.e. corresponding to the Bargmann constraint). How about the 
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case restricted on a subsymplectic manifold in the space B?^ (i.e. corresponding to 
the C. Neumann constraint)? This is a difficult problem. Although r-matrix works 
out and there is no answer about the algebro-geometric solution up to now. 

From section 5, we know that the c-Toda (or r-Toda) flow and the c-CKdV (or r- 
Toda) flow share the same r-matrix as well as the common Lax matrix and involutive 
conserved integrals in the whole space R^^ (or on certain symplectic submanifolds in 
R^^). Therefore a further conjecture is: whether any finite dimensional continuous 
Hamiltonian flow can be associated with a finite dimensional discrete symplectic 
map such that they share a common Lax matrix? If it is right, then the discrete 
integrable systems will be mostly enlarged. 
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